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Abstract 

The present paper provides an overview of results obtained in four recent papers by the 
authors. These papers address the problem of intermittency for the Parabolic Anderson 
Model in a time- dependent random medium, describing the evolution of a "reactant" in 
the presence of a "catalyst". Three examples of catalysts are considered: (1) independent 
simple random walks; (2) symmetric exclusion process; (3) symmetric voter model. The 
focus is on the annealed Lyapunov exponents, i.e., the exponential growth rates of the 
successive moments of the reactant. It turns out that these exponents exhibit an interesting 
dependence on the dimension and on the diffusion constant. 
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1 The Parabolic Anderson Model 



1.1 Motivation 

The Parabolic Anderson Model is the partial differential equation 
d 

— u(x, t) = kAu(x, t) + 7£(x, t)u(x, t), x £ Z d , t > 0, (1.1) 
for the R- valued random field 

u = {u(x,t): x£Z d ,t> 0}, (1.2) 

where k € [0, oo) is the diffusion constant, 7 G [0, 00) is the coupling constant, A is the discrete 
Laplacian, acting on u as 



Au(x,t)= Yl Hv,t)-u(x,t)] (1.3) 
is the Euclidian norm), while 



||y-j;||=l 



i = U(x,t): x G Z d , t > 0} (1.4) 

is an R-valued random field that evolves with time and that drives the equation. As initial 
condition for (jl.ip we take 

n(-,0) = l. (1.5) 

One interpretation of (|1.1|) and (|1.5p comes from population dynamics. Consider a spatially 
homogeneous system of two types of particles, A (catalyst) and B (reactant), subject to: 

(i) yl-particles evolve autonomously, according to a prescribed stationary dynamics given 
by the £-field, with £(x,t) denoting the number of ^4-particles at site x at time t; 

(ii) 5-particles perform independent simple random walks with jump rate 2dn and split into 
two at a rate that is equal to 7 times the number of A-particles present at the same 
location; 

(iii) the initial density of S-particles is 1. 
Then 

u(x,t) = the average number of S-particles at site x at time t . . 

conditioned on the evolution of the ^4-particles. 

It is possible to add that £>-particles die at rate 5 E (0, 00). This amounts to the trivial 
transformation 

u(x,t) u(x,t)e~ st . (1.7) 

What makes (jl.ip particularly interesting is that the two terms in the right-hand side 
compete with each other: the diffusion (of i?-particles) described by kA tends to make u flat, 
while the branching (of S-particles caused by ^4-particles) described by £ tends to make u 
irregular. 



2 



1.2 Intermittency 

We will be interested in the presence or absence of intermittency. Intermittency means that 
for large t the branching dominates, i.e., the li-field develops sparse high peaks in such a way 
that u and its moments are each dominated by their own collection of peaks (see Gartner 
and Konig [10] . Section 1.3). In the quenched situation, i.e., conditional on £, this geometric 
picture of intermittency is well understood for several classes of time-independent random 
potentials £ (see e.g. Sznitman |16] for Poisson clouds and Gartner, Konig and Molchanov [11] 
for i.i.d. potentials with double-exponential and heavier upper tails; Gartner and Konig [10J 
provides an overview). For time- dependent random potentials £, however, such a geometric 
picture is not yet available. Instead one restricts attention to understanding the phenomenon 
of intermittency indirectly by comparing the successive annealed Lyapunov exponents 

A p = lim AJt), p£N, (1.8) 

with 

A p (t) = ilogE([n(0,t)F) 1/p , P £N, t>0, (1.9) 

where E denotes expectation w.r.t. £. One says that the solution u is p-intermittent if 

A p > Ap_i, (1.10) 

and intermittent if (jl.lOp holds for all p £ N \ {1}. 

Carmona and Molchanov [2] succeeded to investigate the annealed Lyapunov exponents, 
and to obtain the qualitative picture of intermittency (in terms of these exponents), for po- 
tentials of the form 

S(x,t) = W x (t), (1-H) 

where {W x (t): x G Z rf , t > 0} denotes a collection of independent Brownian motions. (In 
this case, (jl.ip corresponds to an infinite system of coupled Ito-diffusions.) They showed that 
for d = 1,2 intermittency holds for all k, whereas for d > 3 p- intermittency holds if and 
only if the diffusion constant n is smaller than a critical threshold k p = K p (d, 7) tending to 
infinity as p — > 00. They also studied the asymptotics of the quenched Lyapunov exponent 
in the limit as k j 0, which turns out to be singular. Subsequently, the latter was more 
thoroughly investigated in papers by Carmona, Molchanov and Viens [3j, Carmona, Koralov 
and Molchanov [T], and Cranston, Mountford and Shiga 

In Sections [2HH we consider three different choices for £, namely: 

(1) Independent Simple Random Walks. 

(2) Symmetric Exclusion Process. 

(3) Symmetric Voter Model. 

For each of these examples we study the annealed Lyapunov exponents as a function of d, k 
and 7. Because of their non-Gaussian and non-independent spatial structure, these examples 
require techniques different from those developed for (|1.11|) . Example (1) was studied earlier 
in Kesten and Sidoravicius [12] . We describe their work in Section 12.21 
By the Feynman-Kac formula, the solution of (|1.1|) and (|1.5|) reads 



u(x, t) = E x ( exp 



7 / dsi{X*{s),t-s) 







(1.12) 



3 



where X K is simple random walk on Z with step rate 2dn and E x denotes expectation 
with respect to X K given X K (0) = x. This formula shows that understanding intermittency 
amounts to studying the large deviation behavior of a random walk sampling a time-dependent 
random field. 



2 Independent Simple Random Walks 

In this section we consider the case where £ is a Poisson field of Independent Simple Random 
Walks (ISRW). We first describe the results obtained in Kesten and Sidoravicius [12]. After 
that we describe the refinements of these results obtained in Gartner and den Hollander [B] . 



2.1 Model 

ISRW is the Markov process with state space 

n = (NU{0}f (2.1) 
whose generator acts on cylindrical functions / as 

where the sum runs over oriented bonds between neighboring sites, and 



rj(z) if z^x,y, 

r]{x) — 1 if z = x, (2-3) 

My) + 1 ifz = y> 



i.e., rj xrvy is the configuration obtained from rj by moving a particle from x to y. We choose 
£(• ,0) according to the Poisson product measure with density p E (0, oo), i.e., initially each 
site carries a number of particles that is Poisson distributed with mean p. For this choice, the 
£- field is stationary and reversible in time (see Kipnis and Landim |13j). 

Under ISRW, particles move around independently as simple random walks, stepping at 
rate 1 and choosing from neighboring sites with probability l/2d each. 



2.2 Main theorems 

Kesten and Sidoravicius [12J proved the following. They considered the language of ^4-particles 
and i?-particles from population dynamics, as mentioned in Section \1. 11 and included a death 
rate 5 £ [0, oo) for the 5-particles (recall (jl.7j) ). 

(1) If d = 1, 2, then - for any choice of the parameters - the average number of i?-particles 
per site tends to infinity at a rate that is faster than exponential. 

(2) If d > 3, then - for 7 sufficiently small and 5 sufficiently large - the average number of 
-B-particles per site tends to zero exponentially fast. 

(3) If cZ > 1, then - conditional on the evolution of the ^4-particles - there is a phase 
transition: for small 5 the .B-particles locally survive, while for large 5 they become 
locally extinct. 
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Properties (1) and (2) - which are annealed results - are implied by Theorems 12.21 and 
below, while property (3) - which is a quenched result - is not. The main focus of [T2] is on 
survival versus extinction. The approach in [12], being based on path estimates rather than 
on the Feynman-Kac representation, produces cruder results, but it is more robust against 
variations of the dynamics. 

In Gartner and den Hollander [6] the focus is on the annealed Lyapunov exponents. The- 
orems ETrEH] below are taken from that paper. 

Theorem 2.1. Let d > 1, p, 7 G (0, 00) and p G N. 

(i) For all k G [0,oo), the limit in (jl.8jl exist. 

(ii) If A p (0) < 00, then k — > X p (k) is finite, continuous, non-increasing and convex on [0, 00). 

Let pt(x, y) denote the probability that simple random walk stepping at rate 1 moves from 
x to y in time t. Let 



G,i 



Pt(0,0)dt 



(2.4) 



be the Green function at the origin of simple random walk. 



Theorem 2.2. Let d > 1, p, 7 € (0,oo) and p G N. Then, for all k £ [0, 00), X p (k) < 00 if 
and only if p < l/Gdj- 

It can be shown that if p > l/GdJ, then A p (t) in (jl.9p grows exponentially fast with i, i.e., the 
p-th moment of ii(0, t) grows double exponentially fast with t. The constant in the exponent 
can be computed. 

In the regime p < l/Gd'J, K 1— * A p (k) has the following behavior (see Fig. [T]): 

Theorem 2.3. Let d > 1, p, 7 G (0, 00) and p G N such that p < 1/ G^- 

(i) k 1 — > X p (k) is continuous, strictly decreasing and convex on [0, 00). 

(ii) For k = 0, 

Ap(0)-~' {1/Gd 



(Hi) For k 



00, 



91 (1/G d )- P1 
lim 2dK[X p (K) - fry] = P1 2 G d + l d=3 {2df{ P1 2 pf V 3 



with 



sup 

ll/b=l 



dx\f{x)\ 2 / dy\f(y)\ 



1 



47r||x — y\ 



dx I V/(x) 



(2.5) 

(2.6) 
(2.7) 



d > 4 
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Figure 1: k i— > A p (k) for p = 1,2,3 when p < l/Gd'J for simple random walk in d = 3 and 
d > 4. 

2.3 Discussion 

Theorem 12.21 says that if the catalyst is driven by a recurrent random walk {Gd = oo), then 
it can pile up near the origin and make the reactant grow at an unbounded rate, while if the 
catalyst is driven by a transient random walk {Gd < oo), then small enough moments of the 
reactant grow at a finite rate. We refer to this dichotomy as the strongly catalytic, respectively, 
the weakly catalytic regime. 

Theorem I2.3( i) shows that, even in the weakly catalytic regime, some degree of clumping 
of the catalyst occurs, in that the growth rate of the reactant is > fry, the average medium 
growth rate. As the diffusion constant k of the reactant increases, the effect of the clumping 
of the catalyst on the reactant gradually diminishes, and the growth rate of the reactant 
gradually decreases to fry. 

Theorem I2.3( ii) shows that, again in the weakly catalytic regime, if the reactant stands 
still, then the system is intermittent. Apparently, the successive moments of the reactant are 
sensitive to successive degrees of clumping. By continuity, intermittency persists for small k. 

Theorem I2.3( iii) shows that all Lyapunov exponents decay to fry as K —* oo in the same 
manner when d > 4 but not when d = 3. In fact, in d = 3 intermittency persists for large 
k. It remains open whether the same is true for d > 4. To decide the latter, we need a finer 
asymptotics for d > 4. A large diffusion constant of the reactant hampers localization of the 
reactant around regions where the catalyst clumps, but it is not a priori clear whether this is 
able to destroy intermittency for d > 4. We conjecture: 

Conjecture 2.4. In d = 3, the system is intermittent for all k £ [0, oo). 

Conjecture 2.5. In d > 4, there exists a strictly increasing sequence < k% < K3 < . . . such 
that for p = 2, 3, . . . the system is p -intermittent if and only if k € [0, k p ). 

In words, we conjecture that in d = 3 the curves in Fig. [T] never merge, whereas for d > 4 the 
curves merge successively. 

What is remarkable about the scaling of X p {k) as k — > 00 in (|2.6p is that V3 is the 
variational problem for the so-called polaron model. Here, one considers the quantity 



where a > and (/3(t))t>o is standard Brownian motion on M 3 starting at (3{0) = 0. Donsker 
and Varadhan [5j proved that 



Lieb [13] proved that (|2.7p has a unique maximizer modulo translations and that the centered 
maximizer is radially symmetric, radially non-increasing, strictly positive and smooth. A 
deeper analysis shows that the link between the scaling of \ p (n) for k — > 00 and the scaling 
of the polaron for a — > 00 comes from moderate deviation behavior of £ and large deviation 
behavior of the occupation time measure of X K in (j!.12|) . For details we refer to Gartner and 
den Hollander [B]. 




(2.8) 



lim lim 6{t;a) = 4^7^. 



(2.9) 



6 



3 Symmetric Exclusion Process 



In this section we consider the case where £ is the Symmetric Exclusion Process (SEP) in 
equilibrium. We summarize the results obtained in Gartner, den Hollander and Maillard [7], 

0- 

3.1 Model 

Let p: 7L d x Z d — > [0, 1] be the transition kernel of an irreducible symmetric random walk. 
SEP is the Markov process with state space 

ft = {0,l} zd (3.1) 

whose generator L acts on cylindrical functions / as 

(Lf)( v )= p&y)if(v x ' y )-m], (3.2) 

{x,y}CZ d 



where the sum runs over unoriented bonds between any pair of sites, and 

rf' y {z) -- 



7](z) if z^x,y, 

7](y) if z = x, (3.3) 
jj(x) \iz = y. 



In words, the states of x and y are interchanged along the bond {x,y} at rate p(x,y). We 
choose £(• ,0) according to the Bernoulli product measure with density p £ (0,1). For this 
choice, the £-field is stationary and reversible in time (see Liggett [15J). 

Under SEP, particles move around independently according to the symmetric random walk 
transition kernel p(-,-), but subject to the restriction that no two particles can occupy the 
same site. A special case is simple random walk 

P(x,y) = lf d if f- y|l = 1 ' (3.4) 
I L) otherwise. 

3.2 Main theorems 

Theorem 3.1. Let d > 1, p € (0, 1), 7 G (0, 00) and pgN. 

(i) For all k 6 [0, oo), the limit in (ll.8p exists and is finite. 

(ii) On [0,oo) ; k — > A p (k) is continuous, non-increasing and convex. 

The following dichotomy holds (see Fig. [2]) : 

Theorem 3.2. Let d > I, p £ (0, 1), 7 G (0, 00) and p£N. 

(i) If p(-, •) is recurrent, then X p (k) = 7 for all k £ [0, 00). 

(ii) If p(-,-) is transient, then pj < X p (k) < 7 for all n G [0, 00). Moreover, k \— > X p (k) is 
strictly decreasing with lim^oo X p (k) = pj. Furthermore, p 1— > A p (0) is strictly increasing. 

For transient simple random walk, k t— > X p (k) has the following behavior (similar as in 
Fig. H): 
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Theorem 3.3. Let d>3, p G (0, 1), 7 G (0, 00) and p G N. Assume (f!T4"j) . T/ien 



lim 2dK[X p ( K ) - pj] = p(l - ph 2 G d + l {d=3} (2d) 3 [p(l - p) 7 2 p] 2 P 3 



(3.5) 



mi /i Grf and P3 as defined in (|2.4|) and (|2.7|) . 
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Figure 2: Qualitative picture of k h- ► A p (re) for recurrent, respectively, transient random walk. 
3.3 Discussion 

The intuition behind Theorem 13.21 is the following. If the catalyst is driven by a recurrent 
random walk, then it suffers from "traffic jams", i.e., with not too small a probability there 
is a large region around the origin that the catalyst fully occupies for a long time. Since with 
not too small a probability the simple random walk (driving the reactant) can stay inside this 
large region for the same amount of time, the average growth rate of the reactant at the origin 
is maximal. This phenomenon may be expressed by saying that for recurrent random walk 
clumping of the catalyst dominates the growth of the moments. For transient random walk, on 
the other hand, clumping of the catalyst is present (the growth rate of the reactant is > /ry), 
but it is not dominant (the growth rate of the reactant is < 7). Again, when the reactant 
stands still or moves slowly, the successive moments of the reactant are sensitive to successive 
degrees of clumping of the catalyst. As the diffusion constant k of the reactant increases, the 
effect of the clumping of the catalyst on the reactant gradually diminishes and the growth 
rate of the reactant gradually decreases to pry. 

Theorem 13.31 has the same interpretation as its analogue Theorem I2.3( iii) for ISRW. We 
conjecture that the same behavior occurs for SEP as in Conjectures I2.4H2.5I for ISRW. 

4 Symmetric Voter Model 

In this section we consider the case where £ is the Symmetric Voter Model (SVM) in equi- 
librium, or converging to equilibrium from a product measure. We summarize the results 
obtained in Gartner, den Hollander and Maillard [9]. 

4.1 Model 

As in Section [3l we abbreviate f2 = {0, l} z<i and we let p: 7h d x Z d — » [0, 1] be the transition 
kernel of an irreducible symmetric random walk. The SVM is the Markov process on Q whose 
generator L acts on cylindrical functions / as 




(4.1) 



x,y£Z d 
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where 




In words, site x imposes its state on site y at rate p(x,y). The states and 1 are referred to 
as opinions or, alternatively, as vacancy and particle. Contrary to ISRW and SEP, SVM is a 
non-conservative and non-reversible dynamics: opinions are not preserved. 
We will consider two choices for the starting measure of £: 

{Up, the Bernoulli product measure with density p G (0, 1), ^. 
p p , the equilibrium measure with density p G (0, 1). 

The ergodic properties of the SVM are qualitatively different for recurrent and for transient 
transition kernels. In particular, when •) is recurrent all equilibria are trivial, i.e., \x p = 
(1 — p)5o + p5i, while when p(-, •) is transient there are also non-trivial equilibria, i.e., ergodic 
p p parameterized by the density p. When starting from u p , £(• ,t) converges in law to p p as 
t — > oo. 



4.2 Main theorems 

Theorem 4.1. Let d>l,K& [0, oo), p € (0, 1), 7 G (0, 00) and p G N. 

(i) For all k G [0, 00), the limit in (|1.8|) exists and is finite, and is the same for the two choices 
of starting measure in d^.3[ ) . 

(ii) On k G [0, oo), k — > A p (k) is continuous. 

The following dichotomy holds (see Fig. [3]) : 

Theorem 4.2. Suppose that p(-, ■) has finite variance. Fix p G (0, 1), 7 G (0, 00) and p G N. 
(%) -tf 1 < d < 4, ^/ien A p (k) = 7 /or a// k G [0, 00). 
(ii) If d > 5, i/ien ^7 < X p (k) < 7 /or a// n G [0, 00). 




Figure 3: Qualitative picture of n 1— > A p (k) for symmetric random walk with finite variance in 
d = 1, 2, 3, 4, respectively, d > 5. 



Theorem 4.3. Suppose that p(-,-) has finite variance. Fix p G (0,1) and 7 G (0, 00). // 
d > 5, then p 1— ► A p (0) is strictly increasing. 
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4.3 Discussion 



Theorem 14.21 shows that the Lyapunov exponents exhibit a dichotomy similar to those found 
for ISRW and SEP (see Fig. [3]). The crossover in dimensions is at d = 5 rather than at d = 3. 
Theorem 14.31 shows that the system is intermittent at n = when the Lyapunov exponents 
are nontrivial, which is similar as well. 

We conjecture that the following properties hold, whose analogues for ISRW and SEP are 
known to be true: 

Conjecture 4.4. On [0, oo), k i— >• \ p (k) is strictly decreasing and convex with ]im K — t0O X p (k) = 
fry. 

We close with a conjecture about the scaling behavior for k — ► oo. 
Conjecture 4.5. Let d > 5, p G (0, oo) and p £ N. Assume (13. 4jl . T/ien 

n 2 



lim 2dK[X p (K) - fry] = p(l - p) 7 2 -^ + l {d=5} (2d) 3 



r 5 (4.4) 



with 



POO 

G d = / pt(0,0)dt, 

JO 

/■oo 

/ ip t (0,0)dt, 

JO 



(4.5) 



and 



P 5 = sup 

Il/ll2=l 



dr\f(.rf I d y\f(y)f —^—— - / ,/., j V/(,r ) |" . (1.0) 



5 Concluding remarks 

The theorems listed in Sections [2HH show that the intermittent behavior of the reactant for 
the three types of catalyst exhibits interesting similarities and differences. ISRW, SEP and 
SVM each show a dichotomy of strongly catalytic versus weakly catalytic behavior, for ISRW 
between divergence and convergence of the Lyapunov exponents, for SEP and SVM between 
maximality and non-maximality. Each also shows an interesting dichotomy in the dimension 
for the scaling behavior at large diffusion constants, with d = 3 being critical for ISRW 
and SEP, and d = 5 for SVM. For ISRW and SEP the same polaron term appears in the 
scaling limit, while for SVM an analogous but different polaron-like term appears. Although 
the techniques we use for the three models differ substantially, there is a universal principle 
behind their scaling behavior. See the heuristic explanation offered in [6] and [?]• 

Both ISRW and SEP are conservative and reversible dynamics. The reversibility allows 
for the use of spectral techniques, which play a key role in the analysis. The SVM, on the 
other hand, is a non-conservative and irreversible dynamics. The non-reversibility precludes 
the use of spectral techniques, and this dynamics is therefore considerably harder to handle. 

Both for SEP and SVM, the graphical representation is a powerful tool. For SEP this 
graphical representation builds on random walks, for SVM on coalescing random walks (see 
Liggett [E]). 

The reader is invited to look at the original papers for details. 
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